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Abstract 

This paper is dedicated to provide the global solutions of algebro- 
geometric type for all the equations of a new commuting hierarchy 
containing the Hunter-Saxton (HS) equation. Our main tools include 
the zero curvature method to derive the HS hierarchy, the generalized 
Jacobian variety, the generalized Riemann theta function, the Weyl m- 
fucntions m±(x, t, z), and the pole motion obtained by solving an in- 
verse problem for the Sturm-Liouville equation L(ipi) = —ip" = zyipi. 
Based on these tools and the theory of nonautonomous differential sys- 
tems, topological dynamics and ergodic theory, the algebro-geometric 
solutions are obtained for the entire HS hierarchy 



1 Introduction 

We study here the algebro-geometric solutions for all the equations of a new 
commuting hierarchy containing the Hunter-Saxton (HS) equation, 

Uxxt — UU X xx 2u x Uxxi (!■!) 

where u(x, t) is the function of spatial variable x and time variable t. It arises 
in two different physical contexts in two nonequivalent variational forms 
[TJ [2] . The first is shown to describe the propagation of weakly nonlinear 
orientation waves in a massive nematic liquid crystal director field [1] [2] . The 
second is shown to describe the high frequency limit of the Camassa-Holm 
(CH) equation [TUJQIIGS] 

u t - u xxt + 3uu x = 2u x u xx + uu xxx , (1.2) 
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which was originally introduced in |10t [TT] as model equation for shallow wa- 
ter waves, and obtained independently in |31| with a bi-hamiltonian struc- 
ture. 

The HS equation is a completely integrable system with a bi-hamiltonian 
structure and hence it possesses a Lax pair, an infinite family of commut- 
ing Hamiltonian flows, as well as an associated sequence of conservation 
laws, (Hunter and Zheng [2], Reyes [20]). Traveling waves, asymptotic and 
piecewise smooth solutions of (jl.ip were described in Alber et al El [9] . 
The inverse scattering solutions have been obtained by Beals, Sattinger and 
Szmigielski |19j . Recently, Lenells [23], [23] and also Khesin and Misiolek 
[22] pointed out that it describes the geodesic flow on the homogeneous space 
related to the Virasoro group. Bressan and Constantin [25], also Holden [26] 
constructed a continuous semigroup of weak, dissipative solutions. Yin [27] 
proved the local existence of strong solutions of the periodic HS equation 
and showed that all strong solutions-except space independent solutions- 



blow up in finite time. Gui, Liu and Zhu [28] studied the wave-breaking 




phenomena and global existence. Furthermore, Morozov [29], Sakovich [30] 
and Reyes [20], [2T] investigated (II. ip from a geometric perspective. In [16] . 
we constructed algebro-geometric solutions of the whole HS hierarchy by 
using polynomial recursive and spectral analysis method, due to Gesztesy 



Quasi-periodic solutions (also called algebro-geometric solutions or fi- 
nite gap solutions) of nonlinear equations were originally studied on the 
KdV equation based on the inverse spectral theory and algebro-geometric 
method developed by pioneers such as the authors in Refs.[3]-[3], then in 
[T2| [33] . Roughly speaking, the algebro-geometric approach consists of find- 
ing solutions of HS hierarchy, which are strictly connected to meromorphic 
functions defined on a Riemann surface, in the sense that the zeros and 
the poles of such functions completely determine the solution one is looking 
for. On the other hand, the algebro-geometric setting is directly linked to a 
linear differential operator of Sturm-Liouville type, in particular to its spec- 
trum. In fact, it is well known that if u(x, t) is a solution of HS equation, 
then there exists a positive density function y(x,t) = u xx /2 such that the 
spectrum of the related Sturm-Liouville operator 



does not depend on t [18]. The study of isospectral classes of L leads to the 
determination of the constants of motion for the associated solution u of HS 
hierarchy. 




[12]. 
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Motivated by Zampogni's work p3], in this paper we derive another hi- 
erarchy for HS, linked to the Sturm-Liouville operator L, by using a zero 
curvature method. This approach has a significant difference from our pa- 
per [16]: the algebro-geometric setting of this paper is derives from the 
spectral properties of L, and in particular from the properties of the Weyl 
m-functions. So we give more relevance to the spectral problem than to 
the associated zero curvature relation, i.e. the constants of motion derive 
directly from the isospectral classes of L and not from the zero-curvature 
relation. At the same time, the study of the spectral properties of L is 
carried out using methods of nonautonomous differential equations and the 
classical theory of algebraic curves. 

We pose the problem of finding all the densities y(x) such that Hypothe- 
ses 2.4 hold. It turns out that this problem is equivalent to that of finding 
solutions u(x,t) of the r-th order HS equation having as initial condition a 
solution of the stationary n-th order HS equation. Indeed, there is a surpris- 
ing relation between the Weyl m-function of L and the entries of the matrix 
satisfying the zero-curvature relation, and this fact enables us to find explicit 
algebro-geometric solutions for HS hierarchy, when the initial data lie in an 
isospectral class of L. 

The outline of this paper is as follows. In section 2, we review some 
recent developments in the theory of the inverse Sturm-Liouville problem. 
The material discussed here can be found in \15\ I38j . We use methods of 
nonautonomous differential systems, topological dynamics and ergodic the- 
ory, to characterize all "ergodic potentials p, q, y : R — )■ R" with p, y > 0, 
which constitute a non autonomous differential system (|2.ip . Such a char- 
acterization is carried out by considering the finite poles P\(x), . . . , P n {x) of 
the Weyl m-functions m±(z,x), when interpreted in a "dynamical way". 

In section 3 and section 4, we derive the stationary HS hierarchy and the 
time-dependent HS hierarchy by using a zero curvature method, respectively. 
A solution of the n-th order stationary HS equation is the suitable initial 
condition for solving the r-th order time-dependent HS equation, for every 
fixed n > r G N. 

In section 5, we investigate the relation between the matrix B n of the 
stationary formalism and the Weyl m-functions m± of the Sturm-Liouville 
system. This relation is the key to solving the HS equation. In particular, we 
can see that the entry F n of B n completely determines the Weyl m-functions 
m±, or better, their common meromorphic extension M(P, x) to a Riemann 
surface 1Z. The zeros of F n (when viewed as a polynomial of degree n in z) 
are exactly the poles P%{x) (i = 1, . . . , n) of the function M(P, x). Based on 
this fact and the results of section 2, we give the expression for the solution 
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of the n-th order stationary HS equation in (|5.14p . 

In section 6, we devoted to the time-dependent formalism. We obtain 
the expression for the solution u(x,t r ) of the r-th order HS equation in 
(|6.44p . This solution is of class C°°(IR 2 ), since the poles Pi(x, t) have the 
same regularity properties. Moreover, we study the properties of the Weyl 
m- functions m±. A Riccati type equation with respect to the time variable 
t is obtained. 

In section 7, we investigate the (x, f)-motion of the poles Pi(x, t). It 
turns out that such a motion can be made clear by considering its isomorphic 
image through a generalized Abel map on the generalized Jacobian variety 
of a Riemann surface 1Z of genus n. In particular, we will find that the 
motion on the generalized Jacobian variety can be described by using a basis 
{du\, . . . , duj n ) of holomorphic differentials on 1Z, plus a non-holomorphic 
differential dujQ, in such a way that the x-motion is confined to the non- 
holomorphic coordinate, and defines there a linear function uq(x), while 
the i-motion determines a triangular structure on the first n — 1 of the 
holomorphic coordinates. The image of the t-motion through the generalized 
Abel map is determined by functions Wi(x,t) (i = 1, ... ,n — 1). The last 
holomorphic coordinate co n (x,t) remains implicitly defined in such a way 
that the vector (ojo(x), lo±(x, t), . . . ,co n (x,t)) is contained in a translate To 
of the zero locus of the generalized Riemann theta function 0q(z), for every 
x, t € R. In addition, we give a brief description of how we can build singular 
Riemann surfaces (called "generalized Riemann surfaces" ) , and explain how 
the classical theory can be extended to these new objects. 

In section 8, we show that the elementary symmetric functions of n- 
tuples of distinct points in 1Z can be expressed in terms of a generalized 
Riemann theta function. It turns out that the j-th symmetric function of 
points Pi,...,P n £ 1Z, depends on the partial derivative with respect to 
the coordinate u n -j of the logarithm of wo(Pi(x,t r ), . . . ,P n (x,t r )). As an 
application, we derive the Riemann theta function representation for the 
solution u(x, t) of the HS hierarchy. 

2 Preliminaries 

In this section, we give some important terminology and results about the 
Sturm-Liouville operator, in particular a dichotomy-theoretic approach to 
the inversion problem for such operator. A complete treatment of the results 
in this section can be found in [15j [38] . 

Let A be a compact metric space and {t x ,x £ M} a family of homeo- 
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morphisms of A such that 

(i) To (a) = a, for all a £ A; 

(ii) T x+S {a) = t x o r s (o), for every a £ A and i,s£l; 

(iii) the map r : A x R — )■ j4 : (a, x) i— )• ^(a) is continuous. 

A family {r x } satisfying (i)-(iii) is called a flow on A We fix a {r^j-ergodic 
measure /i on A, and assume that A is the topological support of fi, i.e. 
/J,(V) > for every open set V C A. The triple (^4, {t x }, fi) is called a 
stationary ergodic process. 

Let p, q, y : A — > R be continuous functions with p, y strictly positive; for 
every a £ A, we consider the map (denoted by a(x), with abuse of notation), 



which is equivalent to — (pip'iY + qipi = zyipi, if the map a \-t -^p{t x (o))\ x= q 
is well defined and continuous. 

Let § a {x) be the fundamental matrix solution for (12.ip . 

Definition 2.1 The family of equations (|2.1|) is said to have an exponential 
dichotomy over A, if there are positive constants K, 7 and a continuous 
function a 1— > P a : A — > V = {linear projections P : C 2 — > C 2 } such that the 
following estimates hold: 



It follows from the fact that tr(a(x))=0 for all x £ R that both the im- 
age ImP a and the kernel KerP a are one-dimensional, i.e., can be viewed as 
complex lines in C 2 . The following proposition shows an important charac- 
terization of the spectrum of (|2.ip for /i-a.e. a £ A. A proof can be found 




We study the family of differential equations 




a £ A, z G C, 
(2.1) 



(i) \\^ a (x)P a ^ a (s)- 1 \\ <Ke-^ x ~ s \ x 

(ii) \\$ a {x){I - P a )$ a {s)- l \\ < Ke^ x -^ 



X > s 



X < s. 



in [15] . 
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Proposition 2.2 For ji-a.e. a £ A, the spectrum E a of (|2,ip equals a closed 
set Eel which does not depend on the choice of a. Moreover, 

C\S = {z € C, (|2.1|) has an exponential dichotomy}. 

Let now a £ A and z G C \ S. We define m±(a, z) to be the unique 
complex numbers such that Im P a = Span(l,m+(a, z)) T and Ker P a = 
Span (1, m_(a, z)) T . As usual, we set m + = oo (or m_ = oo) if Im P a = 
Span(0, 1) T (or Ker P a = Span(0, 1) T ). Moreover, it is easily proved that 
Qm+Qz > and Qto-Qz < for every z with Qz ^ 0. For every fixed 
a £ A, the maps z h-> m±(a, z) are analytic in C \ R, and are meromorphic 
in C \ E. They coincide with the classical Weyl m- functions (the definition 
and the properties of the classical Weyl m- functions, see |35tl36|). For every 
z € C, let j3{z) be the upper Lyapunov exponent for (|2.ip . 

Proposition 2.3 Fix a € A. Let I C M be an open interval and suppose 
that (3{z) =0 for a.e. z € /. Then m±(a, z) extend holomorphically through 
1. If h± denote the extensions of m±, then 

)m + (a,z), Qz > 
m_((i, zj, yz < L) 

A proof of this proposition can be found in |15j . 

Throughout all the paper we make the following fundamental assump- 
tions. 

Hypotheses 2.4 (HI) The spectrum S of (|2.ip is a finite union of inter- 
vals, i.e. S = [z 2 , z 3 ] U . . . U [z 2n , oo); 

f-fTS,) The Lyapunov exponent (3{z) vanishes a.e. in the spectrum, i.e. fi(z) = 
/or a.e. z G E. 

Now we give the algebro-geometric structure of this paper. The results 
stated following are proved in |38j . We first look for formula of algebro- 
geometric type for the functions p, q and y. Consider the Riemann surface 
1Z described by the algebraic relation 

2n 

w 2 = z Y[ (z ~ Z m ). 
m=l 

As usual, 1Z is obtained by the union of two Riemann spheres C 2 U {oo} cut 
open along £ and glued together in the usual way, with genus n and exactly 




m_(a, z), Qz > 
m_|_(a, z), Qz < 0. 



6 



2n + 2 ramification points, namely 0,z\, . . . , Z2 n , oo. Let 7r : 7£ — ?■ C U {00} 
be the canonical projection and k(P) be the meromorphic function on 1Z 
defined by 

k(z) = y/(z - Zl) . . . (z ~ Z 2n ), t(P) = Z. 

Then both m+(o, •) and m_(a, •) define a single meromorphic function M a (-) 
on 7£, such that M a = m+ and M a o a = m_, where a is the map inter- 
changing the sheets. M a (P) satisfies the Riccati equation 

Mx + z l / 2 -M 2 =z- 1 / 2 q-z 1 / 2 y, tt(P) = z. 
P 

The function M a has exactly n finite poles P%(x),. . . ,P n (x), one for each 
interval in the resolvent set R\S. For the properties of the function M a and 
a complete discussion about its poles and the algebro-geometric setting, see 

EH. 



If we put T(z) = ~^2( z — 7r (P( 2; )))) then we have 



i=i 

u , p , , p , Q{z) + ^pyk(P) 

M a (P) = m+(a, P) = -j- , tt(P) = z 

l2J (2.2) 
M a ocx(P) = m_(a,i>) = ^-r , 7r(P) = z, 

where Q{z) = ^V" 1 / 2 + q n - X z n -^ 2 + ... + q {x), and Q(tt(P)) = 
y/pyk(Pi), for every i = 1, . . . , n. 

The function M a can be viewed as a function of the variable x by let- 
ting a£ A, and considering the map (x,P) \-t M Tx r a \(P). The same holds 
for p,q,y. It follows that the poles of M T / a \{P) depend on x, hence we 
obtain the "moving poles" P\(x) , . . . , P n (x) . In addition, it is clear that 
Pi(0) = Pi, . . . , P n (0) = P n are the finite poles of the function M a (P). 

Based on this algebro-geometric structure, we are able to find formulas 
for the functions p(x),q(x),y(x) in terms of the finite poles and the zeros of 
M T r a \(P), for every a £ A. For our convenience, we write P%{x) instead of 
7r(Pj(x)), if no confusion arises. 

Theorem 2.5 Suppose that Hypotheses 2.4 hold. Then for every a £ A, 
the finite poles Pj(x) (j = 1, ...,n) of the function M Tx ^{P) satisfy the 
following first-order system of differential equations 



(-irk(P j (x))^P~jx) [m°_(x)-m° + (x)]U7=iPi(x) (l) 
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where m±(x) := m±(T x {a), 0). Moreover, the functions p(x),q(x),y(x) sat- 
isfy the following relations involving the finite poles of the meromorphic func- 
tion M Ta .( a )(P) on TZ, 

2^py T^Pjix)) (-ir +1 [m°(x)-mO(x)]nr=i^(^) 



P k(Pj(x))^P^c) p(x)k(0) 

and 



2n n 2 1 \ 

q(x) = y(x) (E^- 2 E + ^(x) + 



where 



qn {x) = ^ = ^)^l]n([ TO o (a;) _ TO o (a;)]I Jp. (x) 



8=1 



1 V /M^IL - Prix)) p[ V^MII^W " 

The proof of this theorem can be found in [38] . 

What is of great relevance is the fact that this problem can be inverted, 
i.e. the following holds. 

Theorem 2.6 Let z\ < z 2 < . . . < Z2 n be distinct positive real numbers, and 
let TZ be the Riemann surface genus of n described by the following algebraic 
relation 

w 2 = z(z - zi)(z - z 2 ) ■ ■ ■ (z - z 2n ). 

Let P%, . . . ,P n be points on TZ such that z 2 j-i < Pj < z 2 j (1 < j <n). Then 
there exists a stationary ergodic process (A, {t x },/j.) together with functions 
p,q,y : A — >■ M, such that the spectrum of the family 

V2 J V q( T x(a)) - zy{T x (a)) J \ v 2 J 

(2-4) 

of differential equations does not depend on the choice of a £ A and has the 
form 

£ = [Z2,z 3 ] U [zi,z 5 ] U • • • U [z 2n ,oo). (2.5) 

Moreover, (3(z) = 0, for Lebesgue a.e. zeS. The functions p,q,y satisfy 
relations analogous to those of Theorem 2.5. 



S 



The proof of this theorem can be found in [38J . The importance of Theorem 
2.6 lies in the fact that we can generate infinitely many stationary ergodic 
processes and functions p, q, y such that the above conclusions hold. For 
more explanations about this importance, see again [38J. 



3 The stationary HS hierarchy 

In this section, we derive the stationary HS hierarchy and the corresponding 
sequence of zero-curvature pairs by using a polynomial recursion formalism. 
We will use the ergodic-dynamical structure of Section 2 in the following. 

Let C be the set of all the positive uniformly continuous bounded func- 
tions / : R — > R. We equip L with the compact open topology. Let {t x } 
be the Bebutov flow on £, that is, r x (f){-) = f(x + •), for every f £ C. 
Moreover, we define A = Hull(yo) = cls{T x (yo), x G R, yo G £}. One infers 
that A is compact. For every A £ A, let y(A) = A(0). We consider the 
family of equations 

*•-{$)' s)UH** 

(3.1) 

which is equivalent to — ip" = zyip\. 

Now we start with the following 2x2 matrix isospectral problem 

1> x = A(x,z)1> (3.2) 

and an auxiliary problem 

A n = B n (z)il>, (3.3) 
where A(x, z) is a Sturm-Liouville matrix, defined by 

*"■*>-(-«(.) I) {3A) 

and B n {z) is a two dimensional matrix- valued linear differential operator, 
defined by 

*w-(5w r «,w !) ) 26CU0} ' n6N ' (3 ' 5) 

assuming F n , G n and H n to be polynomials of degree at most n with C°° 
coefficients with respect to x. The compatibility condition between f)3.2|) 
and (|3.3p yields the stationary zero-curvature equation 

-B n>x + [A,B n ] = 0, (3.6) 
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namely 



F n , x = 2zG n , (3.7) 

H n>x = -2zyG n , (3.8) 

Gn,x = H n -yF n . (3.9) 

From (|3.7p -(3.9). a direct calculation shows that 

—det(B n (z,x)) = -- — (zG n (z,x) 2 - F n (z,x)H n (z,x)^j =0. (3.10) 

Hence, zG 2 n — F n H n is ^-independent implying 

zG 2 n -F n H n = k 2 (z), (3.11) 

where the integration constant k 2 (z) is a polynomial of degree 2n with re- 
spect to z. If {z m } m= i t ... i 2n denote its zeros, then 

2n 

k 2 (z) = (z ~ Z m ), {2 m } m =l,.,2n £ C. (3.12) 
m=l 

Next, we introduce the following polynomial F n (z),G n (z) and H n (z) 
with respect to the spectral parameter z, 

n 



F n (z) = Y J fiz l , (3-13) 

1=0 
n-l 

G n (z) = ^ 5 ,^, (3.14) 



ffn(«)=X)M I . (3-15) 



z=o 



From (|3.7p . (3.8) and ()3.11|) . we obtain that /o, /to is a constant and /o/io 
— EI z m ■ For our convenience in the following, we choose fo such that 



fo = -h = ~yY\ z m- (3-16) 

Without loss of generality, let fo = ho = 1 to normalize the polynomials F n 
and H n . We keep the same notation for the normalized polynomials, unless 
explicitly stated. 
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We begin our manipulation. 

First, substituting (|3.7p into (3.9) and (3.8), we arrive at 

■yFn,xx ~ zH n + zyF n = 0, (3.17) 

H n , x + yF ntX = 0. (3.18) 
Then, differentiating (|3.17p with respect to x, we obtain 
1 

2 ^n,xxx — %H nx zy x F n zyF nx . (3.19) 
By subtraction and addition of (|3.19p to (|3.18p . we have 

-F n , xxx = -2zyF ntX - zy x F n , (3.20) 

~^F nxxx — 1zH nx zy x F n . (3.21) 

Hence, comparing the coefficients of the same powers in (|3.20p . we have the 
following recursion formalism for the coefficients 

fl,x = -G(4yfl-i,x + (3.22) 

where Q is given by 

I v(x 2 ) dx 2 dx u xeK, v G L°°(R). (3.23) 

-oo J—oo 

It is easy to see that Q is the resolvent of the one-dimensional Laplacian 
operator, that is 

6 =(^) • < 3 - 24 > 



Explicitly, we compute 



/o = l, 

fi = -u + cx, 



1 , (3.25) 

/2 = G{U XX U + -U x ) - UCl + C 2 , 

etc. 

where {q}z<=n C C are integration constants and we have used the assump- 
tion fi(u)\ u=0 = Q,IeI 
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The relation (|3.T|) and (|3.21j) provide the recursion formalism for the 
coefficients g\ and h% of the polynomial G n and H n , one can refer to [16J. 

For fixed n G N, by using (3.9), for I = n we obtain h n = yf n , and 
together with (|3.2ip yields, 

s-HS n (u) = 2u xx f n)X (u) + u xxx f n (u) = 0. (3.26) 

This is the stationary HS equation of order n. 

Hence, we obtain the stationary HS hierarchy. The first equation in the 
hierarchy, obtained by taking n = 1, we compute explicitly, 

s-HSi(u) = — 2u xx u x — u xxx u + u xxx c\ = 0. (3.27) 

It is easy to see that s-HSi(n) = (ci = 0) represents the classical one 
dimensional stationary HS equation. 



4 The time-dependent HS hierarchy 

In this section, we will introduce the time-dependent HS hierarchy. This 
means that u are now considered as functions of both space and time. Let 
r £ N be fixed. We introduce a deformation parameter t r £ M in u, replacing 
u{x) by u(x,t r ), for each equation in the hierarchy. 
Now we consider the Sturm-Liouville matrix 



1 

-zy(x,t r ) 



(4.1) 



and the matrix B r {x, t r ) whose entries are polynomials F r and H r of degree 
r in z, G r of degree r — 1 in z, and with coefficients depending on t r as well. 
Then the compatibility condition yields the zero-curvature equation 

A tr - B r>x + [A, B r ] = 0, r e N, (4.2) 

namely 

-zy tr + H rx + 2zyG r = 0, (4.3) 

F r , x = 2zG r , (4.4) 

G r>x = H r -yF r . (4.5) 

Substituting (4.4) into (4.5) and (4.3), we get 

1 . 
-^F r)XX = zH r - zyF r , (4.6) 
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- zy tr + H r>x + yF r>x = 0. (4.7) 

Then as in the stationary case, differentiating (|4.6p with respect to x, and 
by addition and subtraction to (|4.7p . we arrive at 

^r,xxx = z 2 y tr - 2zyF r>x - zy x F r , (4.8) 

"^Fr.xxx — z yt r ~\~ 2zH rx zy x F r . (4-9) 

Hence, the relation (|4.8p gives the recursion formalism for the coefficients 
fr 

fi = -u + ci, 

^f2,xxx = yt r - 2yfi,x - y x fi, ( 4 - 10 ) 

fi, x = -G(4yfi-i, x + tyxfi-x - ZSipytr)- 

From (4.10), we infer that all the higher order coefficients depend in an 
implicit way on t r , since /2 depends on t r . 

Moreover, the coefficients hi of H r are determined by (|4.9H . Since from 
(4.5), we obtain h r = yf r , then together with (|4.9p yields 

HS r (w) = 2y(x,t r )f r)X (x,t r ) + y x (x,t r )f r (x,t r ) = 0, r > 2. (4.11) 

which is the r-th order HS equation. In addition, in the case r = 1, the 
corresponding equation can also be derived from (|4.9j) . that is, 

HSi(n) = u xxtl + 2u xx u x + «u m - n^ci = 0. (4.12) 

It is clear that HSi (u) = (ci = 0) represents the classical HS equation. 

5 The stationary HS formalism 

In this section we focus our attention on the stationary case. By solving 
the inverse Sturm-Liouville problem, we obtain the relation between the 
Weyl to- functions m± and the entries of the matrix B n of the stationary HS 
hierarchy, which will bring to the determination of the solution of the HS 
hierarchy both in stationary and time-dependent cases. 
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Consider the polynomials F n ,G n and H n as defined before in the sta- 
tionary case, then (|3.11|) gives, 



2;i 



zG 2 n - F n H n = k 2 (z) = J] {z - z m ). (5.1) 



m=l 

We introduce the hyperelliptic curve K n , 

2n 

K n : w 2 - zk 2 (z) = w 2 - z (z - z m ) = 0, {^} m= i,... )2n £ C, (5.2) 

m=l 

which is compactified by joining the point P^ at infinity. The complex 
structure on K n is defined in the usual way [12]. Hence, K n becomes a two- 
sheeted hyperelliptic Riemann surface of genus n (possibly with a singular 
affine part) in a standard manner, denoted by 1Z. 

Let fij, Vj : R — > R (J = 1, . . . , n) be continuous functions with 7^ fj,j, 
vi ^ for i ^ j, and such that 

F„(*) = ( ~n )W+lfc(0) - A*i(«)) (5-3) 

and 

# n (z) = ^~ irfc(0) - i/^x)). (5.4) 

Next, we define the function M : 1Z x R — > C, 

MfPxl fe(^) + V^gn(.) -ffn(z) fp ,_ (55) 

Theorem 5.1 For every fixed P G 1Z such that tt(P) = z6C\S, M(P, x) 
satisfies the following Riccati equation 

M X (P, x) + z 1/2 M(P, x) 2 = -z l/2 y. (5.6) 

Proof. Using (|5.5|) . (3.7), (3.9) and (|3.1ip . a direct calculation shows that 
(EH) is true. □ 



The following theorem describes the nature of the x-motion of the zeros 
{Ati(»)}y=l,...,n of F n . 
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Theorem 5.2 Assume that the zeros {fJ-j{x)}j=i,...,n of F n remain distinct, 
then {fJ>j(%)}j=i,...,n satisfy the system of differential equations, 



, M _ 2{-lTk( H {x))^{x)Y^ =l ^{x) 

Moreover, given initial data satisfying Hj(xo) G [z2j-i, z 2 j] (j = 1> ■■■,n), 
then 

Hj{x) 6 [z 2 j-i,z 2 j], j = l,...,n, x£R. (5.8) 
In particular, fJ,j(x) changes sheets whenever it hits z 2 j-i or z 2 j- 

Proof. The derivatives of (|5.3p with respect to x take on 

F n>x ( H ) = (-l)"^ ) ^ JJ (/i . _ ( 5. 9) 

On the other hand, inserting z = /ij into (|3.7|) yields 

FnA^j) = 1^jGn{Hj) = 2 y /JT j k{lJL : j). (5.10) 

Comparing (|5.9p and (|5.10p leads to (|5.7p . To see (|5.8p . one can use the 
local coordinate A = (z — z*) 1 / 2 near a ramification point z* |12(. I16j. □ 



Now we explain the importance of the relation (I5.7P from two aspects. 
First, it shows the behavior of a point ^j(x) when it reaches a ramification 
point of the resolvent interval containing it. In fact, as soon as the point 
Hj(x) reaches a ramification point z*, its derivative becomes zero (k(z*) = 0), 
and it passes from one sheet to the other in the Riemann surface 1Z, i.e, 
/J.j(x) G [z2j-i,Z2j] for every j = l,...,n. Second, from Theorem 2.6, we 
know that for every choice of points Pj £ [z 2 j-\,z 2 j] (j = 1, . . . ,n), there 
exists a meromorphic function P — > M(x, P) on 1Z, having as finite simple 
poles the points Pj{x) (j = 1, . . . , n), satisfying (j2.3j) . which, adapted to the 
present setting (p = 1, q = 0, m°_(x) = —m° + {x) = 1) becomes 



Hence, the points (ij(x) satisfy the system (|5.1ip . as well as the poles of the 
meromorphic function M y on 1Z defined as in (|2.2p . 
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At this point we can use Theorem 2.6 again to prove that there is a 
function y{x) satisfying the assumptions on the spectrum (HI) and (H2), 
and such that the poles Pj(x) of the function M y satisfy n(Pj(x)) = Hj(x), 
for every j = 1, . . . ,n. 

Another important fact from (|5.7p is that fJ>j(x) € C°°(K), which has 
been proved in our latest paper [To] . 

Based on such a y(x), obtained by solving the inverse Sturm-Liouville 
problem, we obtain the following fundamental theorem. 

Theorem 5.3 Let Pj(x) (j = 1, . . • , n) be n distinct points on the Riemann 
surface 1Z satisfying (|5.1ip . and such that Pj(xo) = Pj £ [z2j-i, Z2j] (j = 
1, ... , n). Then there exists a positive continuous function y(x) such that 

1 - 

F n {z) = -_ H(z - ttCP^x))), (5.12) 

v y{x) j =l 

M y (P,x) = M(P,x) = m + {P,x) and M(P, x) o a(P) = m__(P, x), 

(5-13) 

for every x£l and P G TZ. Moreover, with such y(x), Hypotheses 2.4 hold. 

Remark 5.4 The similar formulas describing the nature of the x-motion of 
the points Vj{x) can be found in ]16^ . here we omit the details. In addition, 
the regularity properties, together with the other observations concerning the 
motion which are valid for fJ.j(x), also hold for Vj{x). 

We must emphasize an important difference between our approach in 
[16j and that in the present paper: while in the former approach, we deduce 
the algebro-geometric structure directly from the hierarchy, in the latter 
we start from an assigned algebro-geometric structure and prove that it is 
conserved by the hierarchy. 

Given these preparations, one of the main results of this section, that 
is, the algebro-geometric formula for the solution u(x) of the stationary HS 
hierarchy reads as follows. 

Theorem 5.5 The solution u(x) of the stationary HS hierarchy is 

n Y ^ n 

u' 



m=l Zm 



where the points Pj(x) (j = 1, ... ,n) evolve according to (|5.1ip . 



16 



Proof. Let z = in (|5.5p . we obtain 

m° + (x) = -l, m°_(x) = +l, (5.15) 

where we used (|3,16p and f)5. 13j) . Hence, take p = 1 in Theorem 2.5, we 
arrive at 

(_-\\n+l n 

V^)= { -j^U p ^)- ( 5 - 16 ) 

Next, we consider the coefficient of z in (|5,12p . 

\n-l n n 



/ 1 (x) = ^ T (^n^( a; ))- ( 5 ' 1? ) 



j=l i=l 



Inserting (|5.16p into (|5.17|) . we get 

AW 



^iII*»(E^)- < 5 - 18 > 

\m=l j=l JV ; 



On the other hand, we note that the coefficient fi(x) in (|3.25p is not nor- 
malized, we find 

u(x)-c 1 = — = fi(x). (5.19) 

11™=! Z r< 



■"m 



Hence 

3=1 jK ' 

The constant c\ can be determined from (|3.1ip . by computing the coefficient 
of z, we obtain 

2n 

Cl = - £ -. (5.21) 

m=l 

6 The time-dependent HS formalism 

The basic problem in an algebro-geometric construction of the solutions of 
the time-dependent HS hierarchy is to solve the time-dependent r-th HS 
equation with a stationary solution of the n-th equation as initial data in 
the hierarchy (recall that r < re). 
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We employ the notations B r , F r , G r , H r , f r , g r , h r to stand for the time- 
dependent quantities, and we will keep the usual notations for the stationary 
quantities. 

Summing up, we are seeking a solution u(x,t r ) of the time-dependent 
algebro-geometric initial value problem 



where 



with entries 



s=0 
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HS r (V) = 2y(x,t r )f r%x (x,t r ) + y x (x,t r )f r (x,t r ) = 0, 
u\t r =to,r = n o, 

S-HS n (li ) = U0,xxx/n + 2U0,xxfn,x = 0, (6.2) 



where 

2y(x,t r ) = u xx (x,t r ). (6.3) 
We start from the zero-curvature equations: 

A tr (x,t r ) - B r)X (x,t r ) + [A(x,t r ),B r (x,t r )] = 0, (6.4) 

- B n:X (x,t r ) + [A(x,t r ),B n (x,t r )\ = 0, (6.5) 



i^rW = ^/ s (x,t r )z s , (6.7) 

s=0 
r-1 

G r (z) = Y,9s(x,t r )z s , (6.8) 



JET P (2r) =^2h s (x,t r )z s . (6.9) 

s=0 

Moreover, it is more convenient for us to rewrite the zero-curvature equa- 
tions (|6.4p and (|6.5p as the following forms, 

-zy tr + H r>x + 2zyG r = 0, (6.10) 

F r ,x = 2zG r , (6.11) 

G rtX = H r - yF r , (6.12) 

and 

F„ x = 2zG n , (6.13) 
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H n>x = -2zyG n , (6.14) 

G n , x = H n -yF n . (6.15) 
From (|6.13p - (|6,15p . we may compute 

^-det(B n (z)) = --^zG n (z) 2 - F n {z)H n {z)) = 0, (6.16) 

and meanwhile Lemma 6.3 gives 

^det(B n (z)) = -~J r (zG n (z) 2 -F n (z)H n (z))=0. (6.17) 

Hence, zG n (z) 2 — F n (z)H n (z) is independent of variables both x and t r , 
which implies 

zG n (z) 2 - F n (z)H n (z) = k 2 {z). (6.18) 

This reveals that the fundamental identity ()3.1ip still holds in the time- 
dependent context. Consequently the hyperelliptic curve fC n is still available 
by (|OD. 

Next, we define 

m + (P,x,t r )= HP) tf^ n{z) , n(P) = z (6.19) 

and _ 

m„(P,x,t r ) = - m +ff n{Z \ ,(P) = z. (6.20) 

The properties of the Weyl m-functions m±(P, x,t r ) are summarized as 
follows. 

Lemma 6.1 The Weyl m-functions m±(P,x,t r ) satisfy the following Ric- 
cati equation, 

M x (P,x,t r ) + z 1/2 M(P,x,t r ) 2 = -z 1/2 y(x,t r ). (6.21) 

Moreover, 

m + (P,x,t r )+m-(P,x,t r ) = 2 ^l G ^ z \ ( 6 .22) 

m + (P,x,t r )-m_{P,x,t r ) = ^fl, (6.23) 

m+(P,x,t r )m-(P,x,tr) = ^44- ( 6 - 24 ) 

F n {z) 
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Proof. The proof of (jIHTj) is identical to (pTEJ) . The relations (p^2l) - (IOID 
are an immediate consequence of the definitions of m± and the fundamental 
identity (16J8D . □ 

Lemma 6.2 Assume that (|6.4p and (|6.5p ZioZd. Let P £ K \ {Poo} 

(x, t r ) £ R 2 . T/ien t/ie function M(P, x, t r ) satisfies the following differential 

equations 

M tr (P,x,t r ) = -z~^H r (z) + 2G r (z)M(P,x,t r ) - Z -^F r (z)M(P,x,t r ) 2 , 

(6.25) 

and 

M tr (P,x,t r ) = (-z-^G r (z) + z- 1 F r (z)M(P,x,t r )) x 

i ~ i ~ ~ (6.26) 

= -z-^H r {z) + z-2yF r {z) + z~ 1 (F r (z)M(P,x,t r )) x . 

Proof. Using f)6. 10j) - ()6. 12|) and (|6.2ip . by a straightforward but rather 
lengthy calculation, we infer 

(d x + 2z^M)(M tr + z~*H r - 2G T M + z~^F r M 2 ) = 0. (6.27) 

Hence 

M tr + z _ 3£r r - 2G r M + Wp r M 2 = C exp( - 2 / z^M dx'^j, (6.28) 

where the left-hand side is meromorphic in a neighborhood of Poo , while the 
right-hand side is meromorphic near Poo only if C — 0. This proves (|6.25p . 
Next, by using (|6.1ip and (I6.2ip . we obtain 

z~%yF r + z~ l (F r M) x = z~3yF r + z~ 1 F r>x M + z- 1 F r M x 

= 2G r M - z~^F r M 2 . (6.29) 

Combining this result with (|6.25p . we conclude that (|6.26p holds. Alterna- 
tively, more efficiently method to prove (|6.25p and (|6,26p can be found in 
our paper [16]. □ 



Next, we study the time evolution of F n , G n and H n by using zero- 
curvature equations (O0]l - ([6TT2l) and ([6TT3]) - (l6TT5l) . 



Lemma 6.3 Assume that fj6.4|) and ()6.5[) hold. Then 

F n ,t r = 2(G n F r - G r F n ), (6.30) 
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zG n> t r — H n F r — H r F n , (6.31) 
H n ,t r = 2(H n G r — G n H r ). (6.32) 
Equations (|6.30p — (|6.32p imply 

-B nA + [B r ,B n }=0. (6.33) 

Proof. Differentiating both sides of (|6,23p with respect to t r leads to 

(m+ - mJ) u = -2k{P)F n , tr F~ 2 . (6.34) 

On the other hand, by (jQ2|) . (^23]) and (^251) . the left-hand side of (^341) 
equals to 

m +: t r — m_j r = 2G r (m + — m_) — z~2 F r (m 2 + — m 2 _) 

= Ak{P)(G r F n -F r G n )F- 2 . (6.35) 

Combining (|6.34p with (|6.35p yields ()6.30p . Similarly, Differentiating both 
sides of (|6.22p with respect to t r gives 

(m+ + m.) tr = 2z\{G nU F n - G n F nU )F~ 2 , (6.36) 

Meanwhile, by (1Q2]) . (^231 and (16^5]) . the left-hand side of (lQ6l equals 
to 

m +jtr + m_ jtr = 2G r (m + + m_) - z~^F r {m\ + m?.) - 2z~^H r 

= -2z^G n F~ 2 F nU + 2z-^F~ l (-H r F n + F r H n $.37) 

Thus, (I6.3ip clearly follows by (|6.36p and (I6.37|) . Next, differentiating both 
sides of (|6.24p with respect to t r yields 

(m+m_) tr = {H nU F n - H n F n)tr )F- 2 . (6.38) 

By using (lfT22]) . (lfT24"]l and flfT25]) . we compute the left-hand side of (jOHjl . 
obtaining 

(m + m_)t r = — z~2 H r {rn + + m_) + 4G r m + m_ 
_i ~ 

— z 2 F r m+m-(m + + m_) 
= 2(G r H n -H r G n )F- 1 -H n F n , tr F~ 2 , (6.39) 

and hence (|6.32p holds. Finally, a direct calculation shows that (|6.30j) - (|6.32j) 



are equivalent to (I6.33p . □ 

The properties of the x-motion and i r -motion of the poles Pj(x, t r ) now 
reads as follows. 
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Theorem 6.4 Assume that (|6.4p and (j6.5[) /ioW. Then, for every j 
l,...,n, 



and 

Pj,t r {x,t r )) 



2(-l) n k(P,(x, t^FrjPjjx, t r )) U7=l Kfa M 

fe(0) ^ *r) n^(^i(^ *r) - Pi(.X, t r )) 
FrjPjMl (xt) 

Proof. It suffices to focus on (16.4ip . since the proof procedure for (16.40P is 
analogous to (|5.7p . Differentiating on both sides of ()5. 12j) with respect to t r 
yields 

F n ,t r {P j ) = -^X[{Pj-Pi). (6.42) 
On the other hand, considering (|6.18j) . we compute ()6.30p at Pj, 



F n ,t r (Pj) = 2F r (Pj ) G n (Pj ) = 2F r (Pj ) . (6.43) 

V Pj 

Hence, combining (pT4"2|) and (j6lE3|) leads to (|6"H|) . □ 



Remark 6.5 ^4 closer look at Theorem 6.4 reveals that the pole motion (both 
for the x-motion and the t r -motion) can be determined by solving only first 
order differential equations 

it = » rf |* = w». 

where U\ and U2 are bounded continuous functions defined on R. 

Now we shall provide the algebro-geometric formula for time-dependent 
HS solutions u(x,t r ). 

Theorem 6.6 Assume that (]6.4[) and (16. 5p /joW. T/ien i/ie r-t/t order HS 
equation (|4.1ip admits a global solution u(x,t r ) of algebro-geometric type, 
when the initial condition u(x, to) = ^o(^) is given by the solution of the sta- 
tionary HS equation of order n. The function yo(x) = y(x, to) 

— l^xx (■£) *o)/2 
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lies in the isospectral class given by Hypotheses 2.4 for the Sturm- Liouville 
operator Lipi = ipi tXX = -zyipi. 
In particular, 



where the pole motion is completely determined from (|6.4U|) and (|6.4ip . 

Moreover, for every t r £ R, the function yt r (%) = y(x,t r ) = u xx (x,t r )/2 
lies in the same isospectral class as yo(x). That is, time evolution of the 
solutions of the r-th order HS equation define densities lying in the same 
isospectral class, and this isospectral class depends only on the initial data. 

Proof. The proof of ()6.44p is analogous to Theorem 5.5. To show that 
all the solutions define densities lying in the same isospectral class, it is 
sufficient to observe that the ^-evolution of the poles Pj(x,t r ) implies that, 
if we start from a time t r , then the motion t r i— > Pj(x,t r ) remains in the 
resolvent interval \z2j—i, Z2j] (j = 1, Hence, we can apply Theorem 

2.6 to conclude the result. □ 

7 Jacobian flows and pole motion on the general- 
ized Jacobian 

In this section we move our attention to the generalized Jacobian variety 
Jq{JZ) of the Riemann surface 1Z. The aim is to make clearer the structure 
of the solutions u(x,t r ) given in (|6.44p . 

We first give a short description of what a generalized Jacobian is, for 
more details, one can refer to [171 134] . 

Let 1Z denote a hyperelliptic Riemann surface of genus n with a standard 
homology basis a^bi (i = 1 , . . . , n) , and dw\ , . . . , dw n be a basis of normal- 
ized holomorphic differentials on 1Z. The normalized differentials du)i, means 



We construct a singular Riemann surface TZq, by pinching a nonzero ho- 
mology cycle at a ramification point Qq on 1Z. Roughly speaking, a stan- 
dard basis of holomorphic differentials is no longer sufficient to describe the 
structure of the Jacobian variety connecting to this new surface TZq: such a 




(6.44) 
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Jacobian variety is called generalized Jacobian variety, and we denote it by 
To solve inverse problems on generalized Jacobian variety, we need the 

(2) 

normalized differential of the second kind <Iwq q on 1Z having a double pole 
at Qo, with principal part Xq^cIXq , The normalization means 

I dw%> =0, (i = l,...,n). 

Let us give a more precise definition of Jo (ft>) ■ We take the set Symm n+l {JZ) 
of unordered n+l-tuples of points on 7Z, such tuples are called divisors. Two 
divisors T>\ and T>i are identified if T>\ — T>2 is the divisor of a meromorphic 
function on TZq. For a more detailed treatment we refer to [38 . The map Iq 




in+l 



sends each (Pq,P\, . . . ,P n ) G Symm n+1 (JZ) to a complex vector (zq,s) G 
C x C n , where P* is a fixed initial point different from the ramification 
points of 1Z. 

Let Ao C C n+1 be the Z-lattice spanned by all vectors of the form 

(2) 

(dw y Q ' o ,dwi, . . .,dw n ), a = a ,ai, . . . , a n , bi, . . . ,b n . 

Clearly J ao dwq^ = and f ao dwi = (i = 1, . . . , n), where ao is a sufficiently 
small simple closed curve centered at Qq G 7^.0. It turns out that the rank of 
A is 2n. The generalized Jacobian J$(JZ) is defined by, Jq(JV) = C n+1 /A . 
It can be shown that Jq{1Z) = C* x C n . For more details about generalized 
Jacobian, one can see [T7l [38] . 

Given this construction, we can solve the Jacoby inversion problem when 
in presence of nonholomorphic differentials. We define the generalized Rie- 
mann Theta function 

@o(z) = z G(s) + cU6( S ), z = (z , s) G C x C", 

where G(s) is the classical Riemann Theta function associated to the Rie- 
mann surface 1Z of genus n, and uj n is defined in (|8.5p . From the standard 

1 ao is required to bound a disc in 1Z centered at Qo, and the support of ao is required 
to be disjoint from the supports of at and bi (i = 1, . . . , n). 
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theory, we know that every symmetric function of divisors (Po, P±, . . . , P n ) 
can be expressed in terms of a Theta quotient. 

We define the generalized Abel map, Iq : Symm n+1 (1Z) — > Jq(IZ) = 
C" +1 /A , by 

I (P ,P 1 ,...,P n )= ( dw< Q ' dWl ' ■ ■ ■ ' dw ni) mod A 0- 

i=0 

Then it is clear that by setting Po = P*, the generalized Abel map Iq 
defines an isomorphism between the space Symm n (lZ) and a non compact 
subvariety of Jo (TV) , given by the locus of the zeros of the generalized Theta 
function ©o- That is, Iq sends divisors of degree n into vectors in Jq(TZ) of 
the form 

I (P 1 ,...,P n ) = z=(jT \dw%l,d Wl ,. . . , dw n )) G C n+1 

1=1 J p * 

such that @o(z — Aq) = 0, where Ao is the generalized vector of Riemann 
constants. 

The subvariety Io(Symm n (1Z)) C Jq{TZ) is denoted by To- It turns out 
that every symmetric rational function of divisors of degree n can be ex- 
pressed as the restriction of the corresponding function of divisors of degree 
n + 1 by setting Pq = P*, hence it can be written as a theta quotient. We 
will study these facts in more detail in Section 8. 

Now we turn to our case, 1Z is the Riemann surface described by the 
following algebraic relation, 

2n 

w2 ( z ) ~ z II ( z ~ z m) = 0, 

m=l 

(2) 

Qq = 0, and the corresponding nonholomorphic differential is dwQ = dwg^. 
Let Cj = 7r _1 [z2j-i, Z2j] (j = l,...,n), each Cj is a simple closed curve 
in 1Z. The product c\ x C2 x . . . x c n is a real n-torus, and embeds into 
Symm n (lZ), and hence into To through the restricted Abel map. Next, 
we choose Pi G c\ , . . . , P n 6 c n , then the correspondence (Pi, ... , P n ) >->■ 

n ,. Pl 

y / (dwQ, . . . ,dw n -i) maps c\ x . . . x c n into a curvilinear parallelo- 
gram IcC™, hence the remaining coordinate \^ l dw n can be viewed as 
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a nonlinear transcendental function on I of y, / (dwo,dw\, . . . ,dw n -i), 

~1 Jp* 



i=i 



such that y / (dioo, rfioi, . . . , dit%_i, dw n ) € To- 
*=i 

Next, we will study the nature of the motion of the poles Pj(x,t r ) on 
Jo (TV). To this aim, let V = (Pi, . . . , P n ) G C n be an arbitrary n-tuple of 
distinct points. We introduce the following symmetric functions: 

^(V) = (-iyJ2 p ei--- p iv e = (£i,...,£i), l<i<n, (7.1) 

where Aj = {£ e N* | 1 < < . . . < h < n}; 
<T?\P) = (-lf 2~2 Pli'-Pti, £ = (h,...,ii), l<i<n-l, (7.2) 



where = {£ € N» | 1 < *i < . . . < U < n, 4 ^ j}. 

The general form of Lagrange's interpolation theorem then reads as follows. 



Theorem 7.1 (Lagrange interpolation formula) LetV = (Pi, . . . , P n ) 6 
C n 6e a n-tuple of distinct points. Then for every k = 1, . . . ,n + 1, i = 
0, . . . , n — 1 and j = 1, . . . , n, we have 

n pk—1 

E tt L p^ fw = *fc,n-< - w- ( 7 - 3 ) 

The simplest Lagrange interpolation formula reads in the case i = 0, 

n pk—1 

£iwW.> = *'"' 

For use in present paper, we recall some important properties of the 
functions and q. These results can be found in [T2], here we omit the 
proofs. 

Lemma 7.2 Let V = (Pi,...,P ra ) £ C n be a n-tuple of distinct points. 
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Then 



(i) + Pjo-Mp) = a^V), i = 0, . . . , n - 1, j = 1, . . . , n. 

r 

(it) ^2<ir-i(P)Pi = aj?\P), r = 0,...,n, j = l,...,n. 



i=0 
r-l 



i 



z- Pi 



i=0 ~ ■? i=0 

We define the (non normalized) differentials 



v s-2 



-cfe, s = 1, . . . ,n. 



2y/l k(z) 

Note that duo is a differential of the second kind, having a double pole at 
0, with coefficient (VvTT^m) °f corresponding principal part. 
Moreover, let 



n r Pi(x,t r ) 



i=l 



doj s -±, s = 1, 



, n. 



p* 



Then differentiating each uj s -\ with respect to x, we have 

8lo s -i(x, t r ) _ ^ Pj 2 { x ^r) dPj(x,t r 



Ox 



^2^/P j (s,t r )A;(P i (a;,t r )) 



(7.4) 



Substituting (|6.40p into (|7.4p and taking account into Theorem 7.1 yields 



-1 



^ fc(0) ^ E[^i( P i(^>*r) - PiOMr)) fc(0) 

which implies the following relation: 



5 s ,i, (7.5) 



dx 



1 



Mo) 



1°. 



Hence we obtain 

u s -i(x,t r ) 



c (t r ) 

^ C-s— 1 Xpr) i 



fc(0) 



, s = 1, 

s = 2, . . . , n. 



x, s = 1, 



(7.6) 



(7.7) 



s = 2, 



, n. 
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This result shows that the x-motion is constant with respect to n — 1 co- 
ordinates, while it is linear with respect to x on the remaining coordinate, 
which corresponds to the differential of the second kind. 

Now we investigate the i r -motion of the poles Pj(x,t r ). Recall ([5. 12|) : 

1 n 

F n (z) = r == l[(z-ir(P i (x,U))). 

By the construction of F r (z), we know that for every r = l,...,n, the 
polynomial F r (z) can be obtained by truncating the polynomial F n (z) at 
the degree r, that is 



F r (z) = F n (z) - z 



r+l 



F n (z) 
7 r+l 



(7.8) 



where [ ] p denotes the polynomial part. Hence, we conclude that 

P[ +1 (*,* 



F r (P j (x,t r )) = -Pj +1 (x,t 



Fn(Pj {%■> ^r)) 



Inserting (|7.9p into (|6.4ip yields 



Pj,t r {x,t r ) 



PUx,t r ) 



y/y{x,t r ) 



(7.9) 



y/y{x,t r ) 

PJ(x,t r )ai- ? 2 r _ 1 (P)fc(P ; -(x,t T .))VP J ( a ;,t r ) 
Yl(Pj(x,t r ) -Pi(x,t r )) 



(7.10) 



As before, differentiating each u s -i with respect to t r gives 

doJ s ^( X ,t r ) Pj'\x,t r ) dPAx,t r ) 



dt r 



f^ 1 2^/PJx~t^k(P j (x,t r )) dt r ' 



(7.11) 



and insertion of (|7.10p into (|7.1ip yields the beautiful relation 



du s -i(x,t r ) 
dt r 



A ^ +r - 2 (x,t r )agl r - 1 (^) 

pif[(Pj(x,tr) ~ Pi(x,tr)y 



(7.12) 



where s = 1, . . . , n and r = 1, . . . , n — 1. 
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After a short computation, we arrive at 

du a -i = 5 a -i tl dtr " 2^ 2^ f=FT- — — dt r 

1=0 \j=l H^i ~*V J 



(7.13) 



Next, we introduce the notation a r = n — r + 1. Then f|T. 13|) gives: 
For every s = 0, . . . , a r — 1, 



For jfe = 0, . . . , r - 2, 

If we set 
and 

then we obtain 



doj s = 5 S idt r . 



n— 1 / n rtk+h. 



ST I 5 ^n-r-h-ljV) \ 



n pk+h 

n h+k (v) = ff — — ^ 



re-1 



h=n—r 

As long as k + /i does not reach the value n — 1, then 'Hh+ki'P) = 0. Hence, 
for every fixed r, we have the following triangular structure for the pole 
motion: 

( 1 



dco s (x, t r 



fe(0) 



0, 



s = 0, 
a = l, 

5 2^ • • • > Cfcy 1^ 



(7.14) 



n-1 



^ H h+ k{V)dt r>h , s = a r + k, fc = 0,...,r-2. 



h=n—r 



A closer look at (|7.14p implies the main differences between our formulas 
for the pole motion in the generalized Jacobian and those in [7]: 
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1. In our context, time motion is confined to the holomorphic coordinates, 
i.e. 0Ji(x,t r ) = wi(t r ), . . . , oj n -i(x,t r ) = co n -i(t r ), while the x-motion 
evolves only in the meromorphic one, namely u)o(x,t r ) = ojq(x); this 
shows a complete separation between spatial and time motions. 

2. For every r = 1, . . . , n — 1, the motion is linear with respect to x, with 
no need of any linearizing change of variables. 

3. The classical one dimensional HS equation corresponds to the case 
r = 1, then the motion on the generalized Jacobian is remarkably 
simple: 

1 

-x + xo, s = 0, 



u) s (x,ti) 



fc(0) 

h+Xi, s = 1 ; 

0, otherwise, 

where xo and Xi are constant phases. 



8 Theta function representation for the solution 
u(x, t r ) on the generalized Jacobian 

In our final section we present expressions for all the elementary symmetric 
functions of n distinct points on 1Z in terms of Riemann theta function. 
In particular, we obtain the theta function representation for the solution 
u(x,t r ) of the HS hierarchy. 

Let 2 

dUk - 1 = 2^W) dZl k = l,...,n. (8.1) 

We note that for k > 1, dojk-i is a non normalized holomorphic differential 
on 1Z. While k = 1, dujQ is a differential of the second kind, having a 
double pole at with principal part (1/ \/0 Zm)^o d\o, in terms of the 
local coordinate Ao = z 1 / 2 near 0. 

Moreover, we denote by dw\, . . . , dw g a normalized basis of holomorphic 
differentials on 1Z, and by dwQ the normalized differential of the second kind 
having a double pole at with principal part Aq 2 d\o, in terms of the local 
coordinate Ao = z 1 ! 2 near 0. 

Let V = (Pi, . . . , P n ), where Pi (i = 1, . . . , n) are distinct points on 1Z. 
For our convenience, we write P{ instead of 7r(Pj) if no confusion arises. 
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Define the variables a%, . . . , a n , such that 
dPj _ 2Pfa { il k {V)k{ P] ) 



j,k = l,...,n, (8.2) 



where Pj = Pj (a\, . . . , a n ). Next, we give the explicit formulas for a±, . . . ,a r 
Combining (|8.ip and (|8.2p . we obtain 

d d AA rPj(ai,.,<x n ) z k-2 



A Pj' 2 9Pj = " Pj-^nUV) 



= £fc,fc — Sn-fe+l^fe,n+l = 1- (8.3) 

Hence, we have 

a fc -a° = V/ du k -x = fe = l,...,n. (8.4) 



V] / ^fc_ 

.7=1 ^ 



Next, we introduce the additional function 

™ /•Pj(ai,...,a n ) z >i-l 

^ 1 = gi, vsfc£) dz = Wn(n (8 ' 5) 

Differentiating (|8.5p with respect to on both sides gives 

dot n P™ 

= Yl ft (p] _p i ) <7 n-fc( 7? ) = <Wl,fc - = -^-fe+l(P). 

(8.6) 

Equation (18.61) shows that all the symmetric functions of the poles P\ , . . . , P n , 
which are restrictions of symmetric functions of n + 1 points in Jq(1Z), can 
be determined by differentiating the n-th non normalized Abel coordinate 
with respect to the {k — l)-th one. 

There is a normalizing matrix D £ Sl(n,C) such that 

a = Dw, (8.7) 

where to and w denote the column vectors of the Abel coordinates (wi, . . . , u) n ) 
and (wi, . . . , w n ) respectively. Moreover, there are constants rj\, . . . ,rj n G C, 
such that 



dwo = k(0)duo + Vidwi, 



i=i 



31 



where 

rji = -k(0) / du . 

J 04 

For s,k = l,...,n, a direct calculation yields 



We denote by D = (y rs ) and D~ l = (/3 rs ). Then from we infer 
n-l 



(8.9) 



<9w 



(8.10) 

where /3 s o = 0. 

Now we are intend to determine the dependence of 7: — with respect to 

the classical Riemann theta function &(s), where s = s(V) = I (Pi, ■ ■ ■ , P n ) 
and / : Symm n (lZ) — > J(TZ) denotes the standard Abel map. 
Since w (Pi, . . . , P n ) = -<? Wn ln @(s(V)) (see E]), then combining ((H3J) 
and (EHJ yields 



s=l 



n— 1 



n-l 



y^7n S ^ s + — -d U] JnQ(s(V)) - k(0)ai - V] 77*10* 



.11) 

Then differentiating (|8.1ip with respect to a k , after some computations, we 
arrive at 



^2f3 sn r] s s n -k+i(V) = ^2 /3 s ,k-lVs + k(0)S 1>k + 



2 



\s=l 



da k uj, 



-ln0(a(P)). 



112) 



Introducing the notation 



£k= ^YlPskVsj , k = l,...,n. 

Next, we are intend to make the clear the meaning of the constants For 
this purpose, consider the (n + 1) X (n + 1) matrix 

/ l/fc(0) -I7l/fc(0) . . . -TJn/HO) \ 





V 



/ 
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It is clear that Dq is a normalizing matrix in the sense that 



\ u n J 



D 



( w \ 



By a short computation, one can infer that the first row of the inverse matrix 
Dq 1 is the vector (fc(0),£i, . . . , £ n ). This implies that the constants are 
those complex numbers such that 



dwo 



i.13) 



fc=i 



Based on the above analysis, we have the following result. 



Lemma 8.1 Let V = (Pi, . . . , P n ) be a n-tuple of distinct points oflZ. Then 
for every k = 1, . . . , n, 



1 



Sn-k+lCP) = T ( 6-1 + HWl,k + 



i) 2 



-lne(s{V)) 



.14) 



where s{V) = I (Pi, ■ ■ ■ , P n ), I ■ Symm n (1Z) — > J(TZ) is the standard Abel 
map, and the constants satisfy (|8.13p . with = 0, at = ^h-X- 

For every n-tuple V = (Pi, . . . , P n ) of distinct points of 1Z, the symmetric 

n 

rational function g : Symm n (lZ) — >• C : V i— )■ JJvr(Pj) is the restriction to 

i=i 

Symm n (lZ), obtained by taking Pq = P* , of a symmetric rational function 
p : Symm n+1 (JZ) — > C. Following [37], every symmetric rational function 
p on Symm n+1 (lZ) defines a meromorphic function on Jq{JZ). Hence, from 
(I8.14p . after some manipulations, we obtain the formula for the symmetric 
function 



7" 



%[h(Pi,...,P n 



(w Q ,wi, . . . ,w n ) - A 



i=i 



u 



Io(Pi, . . . , P n 



{W ,Wl, 



,W r , 



Af 



(8-15) 

where 7 is a constant depending only on the choice of the base point P* and 
the genus n of 1Z, and Ao is the generalized vector of Riemann constants. For 
notational simplicity, we denote the right-hand side of (18.150 by (— l) n G>o('P). 
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We note that (j6.44j) can be rewritten as 



where V(x, t r ) = (Pi(x, t r ), . . . , P n (x, t r )). 

Our main result, the theta function representation of the algebro-geometric 
solution u(x,t r ) for the HS hierarchy now follows from the material prepared 
above. 

Theorem 8.2 The solution u(x, t r ) = u(P\(x, t r ), . . . , P n (x, t r )) of the r-th 
order HS equation, can be written as the following form 

(-l) n+1 ( d 2 ,_ a ,_ f „,_ < w , A ^ 



( fr \ i 

ln@(s(V(x, t T ))) + a - T - ( 8 - 17 ) 

\dUJiUJn J f-f Zi 

1=1 



where V(x,t r ) = ( Pi(x, t r ), . . .,P n (x,t r )), s(V) = I(P x {x,t r ), . . . ,P n (x,t r )), 
and £ n satisfy (|8.13p . and I : Symm n (lZ) — > J(TZ) denotes the standard 
Abel map. 

Proof. The expression (|8.17p is an immediate consequence of Lemma 8.1, 
<|Q)]> and dSHED- □ 

Remark 8.3 Theorem 8.2 shows that the expression (I8.17P depends on the 
theta quotient ()8. 15j) and the partial derivative with respect to the coordinate 
wi of the logarithm of wo(Pi(x,t r ), . . . , P n (x,t r )). Hence, u(x,t r ) can be 
viewed as the restriction of a function, defined on Jq(1Z), to the subvariety 
Tq given by the locus of the zeros of the generalized theta function Qq(z). 



9 Conclusions 

In this paper, we obtained global solutions of algebro-geometric type for all 
the equations of a new commuting hierarchy containing the Hunter-Saxton 
equation. As a main tool we used theta function expressions for all the 
symmetric functions of points P\, . . . ,P n G 1Z. Some of these expressions 
are apparently new. 

On the other hand, the Hunter-Saxton equation belongs to a larger fam- 
ily called Dym-type equation in [TJ 13 [5] , 

u X xt + 2u x u xx + uu xxx — 2ku x = 0, k = constant. 
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One of these equations is a member of the Dym hierarchy that has been 
studied by, amongst others, Kruskal [39J, Cao [3D], Hunter and Zheng [2] 
and Alber et al. 

We remark that although our focus in this paper is on the case k = 0, all 
the arguments presented here can be adapted, with no obvious modifications, 
to study the corresponding equation k ^ 0. As it is observed that by 
substituting y = \u xx into y = \u xx — ^k, then (|4.1ip represents the Dym 
hierarchy and (|4.12p will become the Dym-type equation. The analysis from 
Section 5 to Section 8 can extend line by line to the Dym hierarchy. Hence, 
it is trivial to investigate the algebro-geometric solutions of Dym hierarchy 
again. 
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